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ABSTRACT 

The objective of this project was to study the effect 
that small artificial cavities have on natural convection 
from a horizontal surface. Tests were run with water and 
Freon 113. Data for heat flux as a function of bulk temper- 
ature difference were carefully obtained. These data yielded 
Nusselt number as a function of Grashof number or Rayleigh 
number. All of these were then compared with the data obtained 
by O'Connor [Ref. 3]. Experimental results are presented 
for the heat transfer from horizontal circular disks, with 
and without artificial cavities. 

The data from a mirror finished surface was compared to 
that of previous investigations and the agreement was good. 

The artificial cavities were found to affect the natural 
convection heat transfer from a horizontal surface: the 

heat transfer coefficients in both cases (water and Freon 113) 
were decreased. 
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DESCRIPTION 



Heater Surface Diameter 

Acceleration of Gravity 

Convection Heat Transfer 
Coefficient 

Height of Cylindrical Tank 

Thermal Conductivity of 304 
Stainless Steel 

Thermal Conductivity of Test- 
Fluid Distilled Water or 
Freon 113 

Radius of Cylindrical Tank 

Bulk Temperature of Fluid 

Test Section Surface Temperature 



UNITS 



Inches or ft 



Ft/ 2 
' sec 



BTU/Hr/Ft 2 /°F 
Inches or ft 



BTU/Hr/Ft /°F 



BTU/Hr/Ft/°F 
Inches or ft 



Temperature Difference Between Test 

Surface and Fluid (T - T, ) °F 

s D 

Temperature Measured in Stainless 
Steel Test Section at a Location 
0.227 inch from Test Surface F 

Temperature Measured in Stainless 
Steel Test Section 0.750 inch from 
Test Surface F 



Temperature Difference in Test 

Section (T., - T ) 

1 o 

T + T o 

Film Temperature s b F 

2 

Distance in Stainless Steel 
Section Between Thermocouple 

No. 1 and Thermocouple No. 0 Inches or ft 
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U Subscrint Uncerta inty in Variables Denoted Units of 

p by Subscript, for Example: W T variables denoted 
in T ss by subscript 

S o 

a Fluid Thermal Diffusivity Ft^/sec 



3 Fluid Coefficient of Thermal 

Expansion 

v Fluid Kinematic Viscosity 



1/°F 
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DIMENSIONLESS NUMBERS 



Nu 


Nusselt Number Based on 
Section Diameter 


Test 


h °/ Kf 
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Fluid Praudtl Number 




V /a 


Gr 


Grashof Number Based on 
Surface Diameter 
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g3D 3 6T 

2 


Gr H 


Grashof Number Based on 
Cylindrical Tank Radius 
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Rayleigh Number 




Gr *Pr 
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Power into Heater 




Watt 
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INTRODUCTION 



I . 



A . BACKGROUND 

Some recent experimental data with liquid nitrogen 
showed that drilled cavities in a horizontal disk can improve 
the heat transfer coefficient during free convection [Ref. 1] 
These cavities were approximately 0.015 inches in diameter 
and 0.025 inches deep, and were placed in a copper mirror 
finished, surface. Tests were run with none, one, seven 
and thirteen cavities and the heat transfer coefficient 
continued to improve as the number of cavities were increased 

Free convection heat transfer occurs when a body is 
placed in a fluid at a higher or lower temperature than that 
of the body. The difference in temperature implies that heat 
flows between the fluid and the body and also causes a change 
in the density of the fluid layers in the vicinity of the 
surface. On account of thermal expansion, the fluid at the 
bottom will be lighter than the fluid at the top. Because of 
this there will be a tendency on the part of the fluid to 
redistribute itself, however, this natural tendency on the 
part of the fluid will be inhibited by its own viscosity. 

The difference in density leads to downward flow of the 
heavier fluid and upward flow of the lighter. In both free 
and forced convection, the currents transfer internal energy 
stored in the fluid, however, the intensity of the mixing 
motion is generally greater in forced convection so the heat 
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transfer coefficients are greater than in free convection. 
The heat transfer coefficients in free convection are 
relatively low but good enough for cooling many devices in 
the electrical engineering field: transformers, rectifiers, 

transmission line, electrically heated wires as in incan- 
descent lamps and the heating elements of an electric 
furnace. Free convection heat transfer plays an important 
role in steam radiators, wall of a building, especially on 
the stationary human body in a quiescent atmosphere. A 
knowledge of free convection heat transfer coefficients is 
necessary to understand and determine the heat load on air- 
conditioning or ref regeration equipments. Free convection 
causes the heat losses from pipes or other heated fluids, 
to cool the surfaces of bodies in which heat is generated by 
fission [Ref. 2] in nuclear power applications. 

This experiment examines the effects of surface cavities 
on heat transfer from a horizontal surface with water and 
Freon 113. 

A geometrical configuration should be a first step in 
studying the natural convection in enclosures. The geometry 
and power input situations vary throughout the literature 
of studying the natural convection in enclosures. To start 
this experiment, a circular test section in a cylindrical 
container, with the same size for radius of the tank and 
height of the tank was used. Heat transfer data could be 
compared with existing results, especially with O'Connor's 
[Ref. 3]. 
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The heat transfer across plane vertical and horizontal 
layers of gas, as well as across annular spaces, is most 
interesting from the point of view of applications. In this 
kind of study, we find the problem of insulation by air 
layers, as well as the heat transfer through porous materials 
The boundary layer thickness are no longer small in compar- 
ison with the dimensions of the air space so that the develop 
ment of the flow is often strongly influenced by the shape 
of the boundaries. An unstable stratification is formed if 
the heat flow is directed upward and the lower density fluid 
is located below the cooler fluid (density is higher). This 
unstable stratification must breakdown at a certain value of 
Rayleigh number (Gr>Pr) above which a convective motion must 
be generated [Ref. 4], Jeffreys [Ref. 5] and Low [Ref. 6] 
calculated this limiting value of Rayleigh number and the 
value was 1700 in case of the air layer is bounded on both 
sides by solid walls and a value of 1108 was found with a 
free upper surface. Schmidt and Saunders [Ref. 7] carried 
out an optical investigation in a trough filled with water 
and have confirmed the first of these values very satisfac- 
torily. When the Rayleigh number is above the critical 
value, the convective motion in the fluid takes the form of 
quite regular hexagonals cells (also called Benard cells), 
the fluid moves upward in the center and returns near the 
sides. In case of a horizontal layer with circular side 
walls, a concentric circular pattern of flow was obtained by 
Soberman [Ref. 8]. One method to describe the convective 
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heat transfer across a fluid layer is by an effective 
thermal conductivity Ke , which for a plane layer is defined 
by the equation (Ref. 9]. 

<■ - ¥ < t i - v 

with Ke = thermal conductivity of a solid body which would 
admit the same heat flux as that occurring as a result of 
convection through the fluid layer of thickness 6 under the 
action of a temperature difference (T^ r- Tg) (T^ and T ^ are 
the temperatures of the boundaries of the layer). Let 
K = the true thermal conductivity of the fluid then the ratio 
Ke/K describes the ralative increase in the heat transfer 
because of convection. Ke/K also represent the Nusselt 
number defined with respect to the length 6. The heat transfer 
across horizontal layers of air was measured by Mull and 
Reiher [Ref. 10] and the same measurements performed by 
Schmidt and Saunders [Ref. 7] with water (Pr = 7.1) at a 
Ra = 45000. The transition was quite sharp and could be 
easily determined. The curves for air and water do not 
coincide in a plot of Ke/K versus the Rayleigh number thus 
the true creeping motion can not be assumed, consequently Ra 
does not represent an appropriate dimensionless parameter. 

If a correlation of Nu = C (Gr-Pr m ) n is applied to this 
case one obtains an approximate correlation of the data for 
air and water when m is put equal to 1.65. In the turbulent 
regime, the following equation is obtained 

„ , 1 r , 1/3 

Nu = e/ K = 0.073 (Gr-Pr '” ) 
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1 65 

valid for Gr«Pr 1 > 160,000. In the laminar region the 

scatter of the data is too great to permit a common correla- 
tion. In the case of air the onset of turbulence appears 
less well defined. 

Torrance, et al [Ref. 11] designed an experimental 
model to study laminar flows induced by a simulated fire. 

They used air in a cylindrical enclosure with a heated disc 
in the center. Torrance, et al [Ref. 12] , , compared the 
experimental results with values determined from a numerical 
study. The Nusselt number based on a disc diameter versus 
the Grashof number was plotted. Excellent agreement between 
experimental and theoritical results were obtained. In this 
experimental work, turbulence occurred when the Grashof 

9 

number (based on enclosure height), greater than 1.2 x 10 , 

Lee and Heilman [Ref. 13] have a discussion on fire research 
in their papers. 

The heat transfer from a horizontal circular disc, with 
a mirror finish in a cylindrical enclosure, was recorded by 
Duncan [Ref. 14] who obtained the following relation: 

Nu = 0.31 (Ra) 0,30 

and a comparison of his results with the results obtained 
by Torrance, et al [Ref. 12] theoritical and experimental 
was presented. 

In 1968, P. J. Marto, J. A. Moulson , M. D. Maynard [Ref. 1] 
discussed the nucleate pool boiling of Nitrogen with different 
surface conditions. Pool boiling heat transfer of liquid 
nitrogen from circular, 1-in diameter horizontal disc was 
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studied. The surfaces were: copper and nickel mirror 

finishes, also copper surfaces roughened, grease-coated and 
Teflon-coated. Drilled cavities were cylindrical holes of 
diameter of 0.0043 and 0.015 inch and spark cut conical 
hole of 0.022 inch diameter. Surface roughness and material 
alter the nucleate-boiling curve and the Teflon-coated has 
very little effect while a grease coating significantly 
decreases the nucleate-boiling heat transfer coefficient. 

In this experiment, the effect of artificial cavities on 
both natural-convection and nucleate boiling was determined. 
Many data points were taken in the natural . convection region. 

In 1968, Torrance and Rockett [Ref. 12] made an 

analytical study of the natural convection induced in an 

enclosure by a small list spot centrally located on the floor. 

The enclosure was a circular cylinder, vertically oriented 

with a ratio of 1.0 for height over radius. A Prandtl number 

of 0.7 (air) was assumed and the Grashof number was based on 

the cylindrical height and list spot temperature. The 

equations of fluid flow in axisymmetric cylindrical coordinates 

were solved by numerical method with a computationally stable, 

explicit method. The computation were performed with starting 

conditions from quiescent through the initial transient to 

the fully developed flow. Solutions were obtained for Grashof 

4 4 

number from 4 x 10 to 4 x 10 . The flow patterns and 
temperature fields were illustrated with graphs displaying 
sets of streamlines and isotherms. With an extension of 
Grashof number to 4 X 10^, a periodic vortex shielding 
developed, suggestive of the onset of instability. 
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The cited references and discussion have described what 



has been done in the area of natural convection in enclosures 
B. THESIS OBJECTIVES 



The 


objectives of this work were; 


1 . 


To investigate the effects of drilled cavities on 
the natural convectin heat transfer in water and 
Freon 113. 


2. 


To obtain natural convection heat transfer data from 




a mirror-like surface of circular geometry. The 
results would then be compared to that of Torrance, 
et al [Ref. 12] and also to the correlations of other 
investigators, especially with the results of 
O'Connor [Ref. 3], 



This work would serve as a basis for future studies of 
natural convection heat transfer with different sizes and 
numbers of cavities and with different fluids. 
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II. EXPERIMENTAL APPARATUS 



A. DESIGN CONSIDERATIONS 

The test surface consisted of a stainless steel cylinder 
0.895 inches in diameter installed in the center of a 
cylindrical tank 9.00 inches in diameter and 4.50 inches 
high so that it was flush with the tank floor. There are 
two test surfaces: one blank and one with seven drilled 

cavities. A controllable heat source (using resistance 
heaters) was housed in a copper block, and attached to this 
cylinder. A calibrated resistor in series with the heater 
was used to determine accurately the input power to the 
heater. Eight thermocouples, wired in series were used to 
measure the average bulk temperature of the fluid and one 
removable thermocouple was used to provide a check on this 
temperature. To record the temperature distribution of the 
stainless steel cylinder, four stainless steel sheathed 
thermocouples were equally spaced along the axis of the 
stainless steel test section. The temperatures at these 
locations were obtained using a Dana digital voltmeter. 

The test surface temperature was obtained by linearly 
extrapolating the two nearest thermocouples to the test 
surface (see Fig. 17). 

Figure 1 is a cross sectional view of the tank, test 
section and heater assembly. Figure 2 is an overall view 
of the complete apparatus (including instrumentation). 
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Figure 3 shows a photograph of the two test section cylinders. 
The material for constructing of the tank, test section and 
heater arrangement were available from the apparatus used by 
O'Connor [Ref . 3] . 

The following criteria were used to establish the 
apparatus design: 

(1) The enclosure radius was to be equal to the enclosure 
height. This ratio was convenient to work with and 
the results obtained could be directly compared to 
the experimental and theoritical results obtained 

by Torrance, et al [Ref. 12], Duncan [Ref. 14] and 
O'Connor [Ref . 3] . 

(2) The radiiis of the test surface was to be equal to 
one-tenth of the height of the enclosure, this is 
also the value used in the work of Torrance, et al 
[Ref. 11]. 

(3) The ratio of the depth cavity to diameter was 
approximately two, this corresponds to the value used 
by Moulson [Ref. 15]. 

Insulation was used to minimize the radial heat loss 
due to conduction as much as possible. 

The test section designs used by Mac Kenzie [Ref. 16] 
and Duncan [Ref. 14] were examined. They both used a solid 
cylinder to get one dimentional heat flow. In Duncan's 
work, he used a stainless steel cylinder with a flanged top 
forming the test surface. Duncan [Ref. 14] discussed his 
results and showed the effects of heat transfer due to this 
design. Mac Kenzie [Ref. 16] used a similar cylinder-heater 
assembly in his work. To get one dimentional heat flow in 
cylinder we had to minimize the radial conduction losses. 

The cylinder would not be silver soldered to the heater 
section as O'Connor [Ref. 3] did, but a special design as 
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Figure 1. Cross Sectional View of Tank, Test Section and 
Heater Assembly. 
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Figure 2. Photograph of Complete Apparatus Including Instrumentation. 
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Figure 3. Photograph of the Two Test Sections. 




Figure 4, Blank Test Section. 
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Figure 7. Layout of Seven-Hole Drilled Cavities. 
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shown in (Fig. 6) for a changeable cylinder test section was 
constructed. This test section-heater assembly would be 
rigidly supported at the base of the heater but also it could 
be raised or lowered by screws (see Fig. 6) to keep the test 
surface at the same level with the bottom surface of the tank. 
Because the heat flow through the cylinder would be the basis 
for the calculation of the convection heat-transfer coeffi- 
cient (h), the thermal conductivity of the material used had 
to be accurately known. For this reason and because of its 
ease of machinability , 304 stainless steel was chosen for the 
test surface. 

To minimize the heat loss by conduction through the 
bottom of the tank (since this bottom would be in contact 
with the test section), a material with a very low thermal 
conductivity was chosen. Phenolic was readily available and 
it was used by Mac Kenzie [Ref. 16] in a similar apparatus, 
thus it was known to be strong enough to resist the stresses 
applied in achieving a watertight tank arrangement. The 
stainless steel cylinder was centered in the bottom of the 
tank. To minimize the conduction area, the phenolic was 
machined to a tight fit and beveled. In the phenolic bottom 
plate (Fig. 1) a 0.0625 inch Neoprene rubber "O" ring, 
lubricated with silicone high vacuum grease was installed. 
Plexiglass was used for the tank top and sides. 
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B. TEST SECTION-HEATER ASSEMBLY 



There were two test surfaces used in this experiment: 
one blank and one with seven cylindrical drilled cavities. 

The two-inch stainless steel cylinder (Fig. 3) was machined 
from a one-inch diameter 304 stainless steel rod to a .895 
inch diameter. Four 0.041 inch diameter holes were drilled 
to a depth of 0.375 inch, separated axially approximately 
0.5 inch and 90^* apart for thermocouple installation. The 
stainless steel cylinder was attached to a copper block, by 
means of four screws to facilitate changing test sections 
(see Fig. 6). The copper block served as a housing for the 
heaters. To heat the test surface, five Watlow Firerod 
Cartridge were used, installed in the copper block. A 3.375 
inch section of 2.0 inch diameter stainless steel tubing with 
a 2.0 inch outside and 0.625 inch inside diameter ring 
welded to one end was attached to the copper block and 
heater by two 0.75 inch machine screws, provided as a rigid 
support for the test section-heater assembly. This also 
eliminated the problem of metal to metal contact with the 
test section. In order to minimize the heat loss through 
this support, several holes were cut. 

For seven-hole test section, seven cylindrical holes 
were drilled to a depth of 0.0375 inch and a diameter of 
0.015 inch. Figure 7 is the layout of these seven holes on 
the test surface. 

Several layers of asbestos paper were placed between the 
heater base and stainless steel ring. The stainless steel 
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rings machined to a tight fit to the pipe. These rings were 
bolted to a 12.0 inch diameter phenolic plate, which served 
as a support plate for the entire system. One end of the 
stainless steel pipe was fitted into the stainless steel 
ring. The test section-heater assembly was enclosed in a 
0.125 inch thick steel can with a diameter of 6.25 inch. 

It is important to minimize the heat loss from the test 
section and heater assembly. 

Several layers of asbestos paper were wrapped around the 
heater assembly and a portion of the test section and secured 
with insulating tape. Several layers of Flex-Min K insulation 
were wrapped around the entire assembly, coupletely insulating 
the interior of the stainless steel can. Several layers of 
asbestos paper were put into the hole at the bottom of the 
stainless steel can, between the stainless steel ring holding 
the pipe section and the base plate (phenolic). 

C. INSTRUMENTATION 

1 . Temperature Measurement 

The accuracy of measurement of the temperatures was 
required for this experiment; so, ISA type T copper-constant, 
an thermocouples sheathed in stainless steel (sheath diameter 
diameter = 0.035 inch) were selected. The temperature dis- 
tribution in the stainless steel section was determined using 
four of these thermocouples. The temperature of the surface 
was found using an extrapolating method (see Appendix II, 

Fig. 17). 
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The bulk temperature of the fluid was determined using 
eight Copper Coustantan thermocouples wired in series (four 
pairs) manufactured a Dynatech type of welder and a single 
removable thermocouple to provide a check on this temperature. 
The eight thermocouples were wired in series into an insulated 
switching box so the reading gave an average temperature. 

The sheathed thermocouple leads were securely fastened to the 
bottom phenolic plate of the assembly and the more plexible 
thermocouple extension wire was used to extend the wire into 
the switching box. A reference junction thermocouple was 
inserted in a sealed oil-filled glass tube which was immersed 
in a crushed ice-distilled water bath. This thermocouple was 
also wired into the switching box. 

The output from the thermocouples was read on a Dana 
digital voltmeter accurate to 0.001 milivolt. 

2. Power to Heater 



PRECISION 




Figure 8. Simple Circuit for Measuring Power to Heater 
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In order to measure the power input to the heater, 
a circuit as shown in Figure 8 was constructed. Since R R 
is a function of temperature, we should use a more stable 
resistor (R) (calibrated resistor). With this simple arrange- 
ment, a voltage across heater (V R ) and across resistance 
(V ) were taken. The current through the circuit could be 
determined : 



I H = VR/R (3) 

Knowing R (2.013 ), V„ and V_, , power was calculated 

ii K 

from : 



P = 



Vh 



and from ( 1 ) : 



P = V h (V r /R) (4) 

To get the input power to heater accurately for every 

run, we needed to measure the values of V TJ and V n from (4) 

n K 

only. 
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Ill, 



EXPERIMENTAL PROCEDURE 



A, TEST SURFACE PREPARATION 

1 . Mirror Finish 

A procedure as outlined below was used to prepare 
each test surface; 

The test surface was first sanded with emery paper 
number 1 and 3. After that, it was washed with detergent 
and warm water to remove the previous abrasive also a final 
wash with alcohol. 

In order to get a mirror finish surface the test 
surface was then wet polished on metallurgical wheel covered 
with velvet. It was impregnated with one micron diamond 
dust. The test surface was placed face down on the wheel and 
a low speed of rotation was used. Usually, it took from 
fourty five minutes to an hour to complete the work. During 
the work, alcohol was used to wet the surface at. all time. 
After completion of the polishing, the test surface was washed 
with alcohol again and was dried with a dryer. After work was 
completed, it was wrapped in special paper to prevent con- 
tamination and touching. 

2 . Cylindrical Cavities 

The seven-hole test surface was drilled before the 
mirror finish. The purpose of making a mirror like surface 
to insure that the surface conditions were not responsible 
for the heat-flux of the surface. 
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Cavities were drilled using Sphinx spirec pivot drill, 
in 0.015 inch diameter, The cavity depth was controlled by 
using a fesler gauge on the lathe. A ratio of approximately 
2.5 was desired for cavity depth over diameter. So the 0.015 
inch holes were drilled to a depth between 0.035 and 0.040 
inch . 

The locations of seven holes on the test surface was 
shown in Figure 7: One in the center of the surface and the 

other six equally spaced on a circle whose center was at the 
first hole. The holes were spaced 0.165 inch apart. 

All of the final measurements were taken to determine the 
location of the center of the holes for sheathed thermocouples 
relative to the test surface and the final test diameter 
was measured. An engineering flat and a scale readable up 
to 0.0001 inch were used. 

B. TEST SURFACE INSTALLATION 

Immediately before installation of the seven-hole test 
surface, it was put in an ultrasonic cleaner for ten minutes 
then it was washed with distilled water, rinsed with acetone 
and dried. 

C. APPARATUS ASSEMBLY 

There are two main parts in the apparatus: 1) the 

test sections and heater assembly, 2) the tank. 

The test heater assembly (section III, B) was rigidly 
bolted to the phenolic base. After finishing the insulation 
around the test section and heater assembly, the can was 
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installed with the insulation inside, A phenolic disc used 
as the bottom of the tank was fitted to the test section 
carefully. The level of the test section was adjusted from 
the bottom base of the assembly by screws. Three long 
threaded rods were used to secure the upper plate to the 
phenolic base plate, 120° apart. These three rods also used 
to adjust the test surface in level with with phenolic 
bottom plate. A rubber "0" ring was fitted around the test 
surface and the hole in the bottom plate to prevent leaking. 
For water, rubber gaskets were used at the bottom and top 
of the cylindrical plexiglass to make sure that the enclosure 
did not leak. To secure the top plexiglass plate to the 
phenolic bottom plate of the tank, six threaded rods were 
used, 60^ apart. 

To measure the average bulk temperature of the fluid, 
four small holes were drilled through the top plexiglass of 
the tank to insert four pairs of themocouples all wired in 
series, near the tank wall. One bigger hole was used to 
pour the fluid into or out of the tank, and was also used 
as a hole to insert one removable thermocouple used to provide 
a check on this temperature. It was possible to remove the 
plexiglass top plate and also the cylindrical wall of the 
tank without disturbing the test section in order to change 
the fluid for another run. 

It should be noted that the test surface and phenolic 
bottom plate were cleaned with acetone prior to filling the 
unit. All other surfaces in contact with the fluid were 
cleaned with distilled water. 
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D, FLUID PREPARATION 

The test fluid were distilled water and Freon 113. The 
fluid could be poured into the tank through a one inch 
diameter hole. A siphon was used to fill the tank. 

E , TEST PROCEDURE 

1 . Temperature 

For water, five settings of power input to the heater 
(10, 20, 30, 40, 50 watts) and for Freon 113 seven settings 
(1, 2, 3, 4, 5, 6, 7 watts) were used. Data were taken at 
these settings of power input. At the beginning of each 
test, a lowest power input was set by the method discussed 
in section II, C 2. A long time was allowed to elapse until 
the steady state was reached (usually two to three days for 
each data point). The same procedure was repeated for each 
higher power input setting. When we reached the highest power 
input to the next lower one. When the starting point had 
reached steady state a complete set of data was obtained for 
this run. 

For each setting, we took data for: 1) themocouples 

readings in the stainless steel test section, 2) the bulk 
temperature of the fluid, 3) voltages across heater and 
calibrated resistor. 

Temperature readings were usually taken every hour in the 
daytime only. When the bulk temperature and the temperature 
nearest the heater were changed less than half degree per 
hour, steady state was assumed to have been reached. 
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2 . Determination of Nussolt , Grashof and Rayleigh Numbers 



To determine the Nusselt, Grashof and Rayleigh numbers 
the fluid properties were calculated at the film temperature 

which is the average of the bulk temperature (T, ) and the 

T.+T 
b s 

test surface temperature (T ) i.e. ^ — • The thermal 

conductivity of 304 stainless steel (K ) was calculated 

s s 

using the data from Touloukian [Ref. 17] at the average of 

the temperature of the two nearest tempermocouples to the 

T +T 

test surface i.e, T gs = 0 1 . This corresponds to the average 

2 

temperature of the length Ax = (x^-Xq). 

For each run, the temperature distribution in the stain- 
less steel was found to be linear. For this reason, the 
test surface temperature could be determined by using linear 
extrapolation 



T = B 1 x + B 2 



(5) 



Using the two nearest themocouples readings to the test 
surface, substituting into equation (5) and note that, at 
surface, x = 0, thus: 



T 

surface 



B 



2 



and 



B^ a heat flux at surface 



AT _ (T 1 V 

Ax (x x - x Q ) 



( 6 ) 



(7) 



A simple energy balance at the surface, using Fourier's 
Law of Conduction yields: 
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( 8 ) 



ST 



K 



ss 



ss Sx ^ ^surface ~ ^bulk^ 



The Nusselt number was obtained from the definition: 



and 



Q = 



m hD 

NU = Tr— 

K 



hA(T s - T b ) 



AT 

K A (£±) 
ss kx' 



( 9 ) 



( 10 ) 



or 



h = 



Q 



Q 



A(T S - T b ) A(B 2 - T 4 ) 



From (9), Nu = 



QD 



A(B 2 - T 4 )Kf 



(ID 

( 12 ) 



From (10) and (7), we thus obtained 



Nu = 



B-.DK 

1 ss 



K f < B 2 - T 4> 



(13) 



The Grashof number was calculated using the diameter of 
the heated test surface as the characteristic length: 



Gr = 



g3 D° (T s - T b ) 



(14) 



or 



Gr = D 3 (B - T 4 ) 
v 



(15) 



The Rayleigh number is the product of the Grashof number 
and Praudlt number: 
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Ra = Gr • Pr 



(16) 



A simple calculation for water is contained in Appendix 
B, 1 and for Freon in Appendix B, 2. 

If we use the heigh or radius of the tank as a character- 
istic length then the Grashof number becomes: 

Gr H = Gr R = (I f> h3 (B 2 " T 4 ) (17) 

Since the test section diameter was chosen to be one fifth 
the heigh of the tank thus we have the relation for Grashof 
numbers : 



Gr^ = Gr R = 125 Gr 



(18) 
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TABLE I 



UNCERTAINTIES IN THE EXPERIMENT (For Water) 



Power input 
Settings (Watts) 




Percent Uncertainty 




11,12 


Nusselt 

6 


Number Grashof Number 

6 


Rayleigh Numb* 
6 


20.12 


4 


4 


4 


30.93 


3 


8 


8 


40.93 


3 


2 


3 


50.54 


3 


1 


2 
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TABLE II 



UNCERTAINTIES IN THE EXPERIMENT (Freon, Seven-Hole) 



Power input 
Settings (Watts) 




Percent Uncertainty 




1.16 


Nusselt 

8 


Number Grashof Number 
8 


Rayleigh Numbi 

7 


2.13 


9 


8 


8 


3.13 


7 


9 


9 


4.14 


8 


9 


8 


5.13 


8 


8 


9 


6.10 


8 


7 


9 


7.27 


7 


8 


7 
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TABLE III 



UNCERTAINTIES IN THE EXPERIMENT (Freon, Blank) 



Power input 
Settings (Watts) 




Percent Uncertainty 






Nusselt Number 


Grashof Number Rayleigh Numb< 


1.12 


7 


6 


6 


2.13 


8 


8 


6 


3.13 


7 


6 


6 


4,13 


6 


7 


7 


5.28 


7 


8 


7 


6.51 


7 


7 


6 


7.26 


6 


5 


7 
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IV. RESULTS AND DISCUSSION 



This investigation covers a range of Grashof number (based 

(3 Q 

on test surface diameter) from 1 x 10° to 16 x 10° for water 

with seven-hole test section and for Freon 113 a range from 
Q Q 

4 x 10° to 25 x 10° with seven-hole test section and from 
6 & 

3 x 10 to 23 x 10 with blank test section. A range of 
Prandlt number from 2.58 to 4.65 (with seven-hole test section 
for water and for Freon 113, at a range from 7.65 to 8.84 
(seven-hole test section) and from 7.52 to 8.38 with blank 
test section. A summary of the results can be found from 
Table IV, V and VI. 

Two different surfaces were tested and twenty-nine runs 
were made. The natural convection data points were plotted 
in Figure 9 through 16 and approximated by straight lines. 

In section II, C the bulk temperature of fluid was cal- 
culated. To reach steady state for each input setting, a 
minimum of twenty-four hours was needed. Figure 9 is a 
log-log plot of the Nusselt versus the Grashof number for 
water with seven-hole test surface. For this plot, a vei’y 
good agreement between the heating and cooling phase of 
Nusselt and Grashof number (based on test surface diameter) 
was obtained. 

A brief discussion of the result of water with seven-hole 
test section as compared with the work of O'Connor with 
blank test section [Ref. 3] is needed to have some 
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conclusions on the effect(s) of drilled cavities on heat 
transfer (ref to Figures 9, 10). The log-log plot of the 
Nusselt number vs Grashof number (based on test section 
diameter) was approximated by a straight line which has the 
same slope as that of the least square fit of O'Connor's re- 
sults but with a lower intercept. The work of Mr. O'Connor 
was discussed in the introduction of this experiment. Also, 
a comparison of the result with that of Torrance [Ref. 11] 
is illustrated in Figure 9. A Prandlt number of 0.7 (air) 
was assumed by Torrance and the Grashof number was based on 
cylinder height and list spot temperature. 

For Freon 113 with seven-hole test section, the approxi- 
mate straight line is lower than that of the blank test 
section but the slopes are the same in both cases. A log- 
log plot of the Nusselt number versus the Grashof number for 
seven-hole test section and blank test section is shown in 
Figure 15. 

Some of the data points were taken during the cooling 
phase (when power inputs were decreased). These points were 
denoted by << ’ >> in the Table. For cooling phase at a 
power input setting, to compare with that of the heating 
phase at the same Grashof number the Nusslet number is a 
little bit lower. For Nusselt number versus Rayleigh number, 
the data for cooling phase (correspond to the same power 
input setting) has almost the same Rayleigh number with the 
corresponding point in the upward phase but the Nusselt 
number is lower (for the case of water with seven-hole test 
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section, Figure 9); for Freon 113 the corresponding point 
in the cooling phase has almost the same Ra but the Nu is 
lower. The Nusselt number vs Rayleigh number (log-log plot) 
is plotted in Figure 10 for water with seven-hole test section 
to compare with the result of O'Connor [Ref. 3] . Figure 12 
is a plot for the case of Freon 113. 

The heat transfer correlation which results from this 

experimental data is: 

0 29 

Nu = 0,38 (Ra) ' for water with seven-hole test section 

Nu = 0.48 (Ra)O'^ for Freon with seven-hole test section 

and Nu = 0.53 (Ra)^'^ for Freon with blank test section. 

Torrance [Ref. 12] in his numerical solution, assumed a uniform 

fluid bulk temperature and in his experimental work, used 
an average of the container wall temperature and the floor 
temperature. He obtained a good agreement between his 
experimental and numerical method. In this thesis, an average 
bulk temperature was used (automatically obtained from the 
reading of eight thermocouples inserted in the fluid). 
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SUMMARY OF RESULTS FOR WATER (With 7^Hole Test Section) 
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indicates cooling phase 



SUMMARY OF RESULTS FOR FREON (With 7-Hole Test Section) 
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indicates cooling phase 



SUMMARY OF RESULTS FOR FREON (With Blank Test Section) 
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V, CONCLUSIONS AND RECOMMENDATIONS 



The results lead to the following conclusions: 

1. Cylindrical drilled cavities affect the natural 

convection heat-transfer of water and Freon 113: the 

line of log-long plot of the Nusselt number vs Grashof 
number and Rayleigh number with seven-hole test section 
is lower than that of O'Connor (with blank test section) 
for water, but the slope is almost the same. 

2. The heat transfer line in the log-log plot of the 
Nusselt number versus Grashof number and Rayleigh 
number for Freon 113 is higher than that of water with 
the same test section. 

3. For water with seven-hole test section, a correlation 
which results from this experimental data is: 

0 9Q 0 

Nu = 0,38(Ra) ' for water and Nu = 0.48(Ra) 
for Freon 113 with seven-hole test section and 
Nu = 0 . 53 ( Ra)° ‘ 30 

for Freon 113 with blank test section 
The following recommendations are made for future 
work : 

1, Continue the present studies with different size and 
number of cavities (say one, thirteen, and ninety six 
cavities) and compare the results with this experimental 
work . 

2. Using a power setting of at least eight watts for the 
case of Freon 113 to get a boiling curve and the result 
would be compared to that of Mr. Moulson [Ref. 15]. 

3, Using the same test sections but different fluids to 
investigate the effects of changing fluid. 

4, Study the effects of different cavity geometeries 
(say cylindrical, conical, reentrant). 

5. Find a method to get the bubbles out of the tank if 
the power input is greater than 40 watts in case of 
water to get a good result on heat transfer. 

6. Rework this experimental stady with the flow visual- 
ization technique and compare the results with the 
results of this study. 
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Figure 15 
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Figure 16 
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APPENDIX A 



THERMOCOUPLE CALIBRATION PROCEDURE 

The calibration of the sheathed and bulk temperature 
thermocouples was a major part of this experiment since all 
the results and discussion based on the accuracy of the 
determination of the test surface and fluid temperatures. 

All the wiring and recording instruments would be used in 
the experiment were not changed after the calibration was 
completed . 

A Rosemount type of calibration system was used with a 
constant temperature oil path. The ice junction thermocouple 
(sealed in a glass tube, filled with oil) was placed in a 
stainless steel dewar filled with ice water and circulated 
with a pump. The nine bulk temperature thermocouples and 
the four sheathed thermocouples were suspended about three 
about three inches into the oil bath. The calibration was 
conducted over a range from 80^F to 324°F and a maximum error 
for the thermometer for this temperature range was 0.005°F. 

A Platinum Resistance thermometer with a commutation bridge 
was used as a standard. 

A comparison of the calibration results with standard 
thermocouple tables: for the sheathed thermocouples the 

result was good with deviation from the table values in a 
range from -0.012 mv to +0.005 mv and for the bulk temperature 
in a range from -0.009 mv to +0.001 mv. 
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The above results led us to use the standard tables with 



an uncertainty in the sheathed thermocouples reading of 
0.5°F and 0.5°F in the bulk temperature reading. 
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APPENDIX B 



SAMPLE CALCULATIONS 



I , FOR WATER 

These calculation are for data recorded on 20 December 
1974 for an input of 11.12 watts and with seven-hole test 
surface . 

1 . Determination of the Surface Temperature 

To determine the temperature of the test surface a 
graphical method was used. A plot of the temperature versus 
the thermocouple distance from the test surface was made for 
each run (corresponding to a proper setting of power input), 
as shown in Figure 17 using the two nearest sheathed thermo- 
couples to the test surface. Using a linear extrapolation 
to obtain the surface temperature an error within 0.5°F 
could be obtained. 

2. Calculation of the Nusselt Number 

Using the equation (13) of section III, E2 to calculate 
the Nusselt number: 



Nu 



B l DK ss 

K f (B 2 ~ T 4 ) 



( II-A) 



The value of K was obtained based on an average 
s s 

temperature of the cylinder over the length Ax using the two 
nearest temperatures to the test surface in the stainless 
steel cylinder. The values of Kss were obtained using the 
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Figure 17 



DETERMINATION OF SURFACE TEMPERATURE 



Power Input Setting: 11.12 Watts 

Bulk Temperature: 85.60 °F 
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data from Touloukian [Ref, 17], Fluid properties were 

T T 

evaluated at the film temperature i,e. = — S — 2 



(section III, E2) using the data from the Table VI, Appendix 
C, from Kreith [Ref. 18]. For 11.12 Watts input setting, 
the following were obtained: 

T Q = 120 . 04°F 
T x = 139 . 54°F 
T 4 = 85 . 00°F 

D = 0.07458 ft. 

Ax = 0.04183 ft. 

= 566.14 Btu/hr-ft 3 
B 2 = 110 . 15°F 
Tss = 129 . 79°F 
T f = 97 . 88°F 



Kss = 9.099 Btu/hr-f t-°F 
K f = .363 Btu/hr-f t-°F 
gB/ v 2 = 111.000 x 10 6 1/ f t 3 -or 
Pr = 4.65 



and from (II-A) 



(466. 14) (0.07458) (9. 099) 
WU .363(110.15-85.60) 



35.49 



3. Calculation of the Grashof Number 

Using equation (15) of section III, E2 to calculate the 
Grashof number : 



Gr = (g$/ v 2) D 3 (B 2 - T 4 ) (15) 

The values gS/ v 2 were calculated from Appendix C, from 
Kreith [Ref . IS] . 
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For 11,12 watts input setting, the Grashof number is; 

Gr = 111,00 x 10 6 (0.07458) 3 (110,15 - 85,60) = 1.13 x 10 6 

4 , Calculation of the Rayleigh Number 

Using equation (16) of section III, E2 , the Rayleigh 
number is the product of the Grashof number and Prandlt 
number : 

Ra = Gr • Pr (19) 

using (12) the Rayleigh number is 

R = 1.13 x 10 6 (4 . 65) = 5.25 x 10 6 

II. FOR FREON 113 

The following is a sample calculation of the Nusselt, 
Grashof and Rayleigh numbers for data recorded on 19 March 
1975 for 1.16 watt of power input setting and with seven- 
hole test surface. 

1 , Determination of' the Surface Temperature 

Using the same method in 1-1, the test surface temperature 
was obtained in Figure 18. 

2 . Calculation of the Nusselt Number 

Using the quation (II-l) of section III, E2 to calculate 
the Nusselt number: 

„ (B 1 )(D)(K SS ) 

m K f (B 2 - T 4 ) (II-l) 

The values of K were obtained using the data from 

s s 

Touloukian [Ref, 17], Fluid properties were evaluated at the 

T + T b 

film temperature T^ = K - ^ (section III, E2) using the 
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Figure 18 



DETERMINATION OF SURFACE TEMPERATURE 



Power Input Setting: 1.16 Watt 

Bulk Temperature: 76.17 °F 
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data from Table n and Table 13 from Thermophysical Properties 



of Refrigerant 113 [Ref. 19], For 1.16 watt of power input 
setting, the following were obtained: 

T q = 84 . 23°F 
T x = 86 , 12°F 
T 4 = 76 . 17°F 
D = 0.07458 ft. 

B 1 = 45.18 Btu/hr-f t 2 

B 2 = 83 , 27°F 

Tss = 85 . 17°F 

T f - 79 . 72°F 

Kss = 8.81 Btu/hr-ft-°F 

y _= 1.592 lbm/ft-hr 

K f =0.0432 Btu/hr-f t-°F 

Cp = .229 Btu/lbm-°F 

p = 97.47 lbm/ft 3 

The values of $ were calculated by definition from the 
thermo-dynamic properties [Ref. 20], 



6 



1 dp n _ 1 / p2 - pi . 

p dT ; p p k T 2 - T 1 ; 




( 20 ) 



with 



v f 



v f 2 + V fl 



For 




1,16 watt of power input setting, we have: 

6 = 8.28 x 10" 4 '*'/ <, R using a temperature difference 
- 10°F in equation (20). 
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